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Abstract. In this article we study the time evolution of an interacting field theoretical 
system, i.e. (/>*-field theory in 2+1 space-time dimensions, on the basis of the Kadanoff- 
bJQ, Baym equations for a spatially homogeneous system including the self-consistent tad- 

^ ' pole and sunset self-energies. We find that equilibration is achieved only by inclusion of 

I the sunset self-energy. Simultaneously, the time evolution of the scalar particle spectral 

function is studied for various initial states. We also compare associated solutions of the 
corresponding Boltzmann equation to the full Kadanoff-Baym theory. This comparison 
shows that a consistent inclusion of the spectral function has a significant impact on 
the equilibration rates only if the width of the spectral function becomes larger than 
' 1/3 of the particle mass. Furthermore, based on these findings, the conventional trans- 

, port of particles in the on-shell quasiparticle limit is extended to particles of finite life 

' time by means of a dynamical spectral function A{X,p, M^). The ofi'-shell propagation 

OO I is implemented in the Hadron-String-Dynamics (HSD) transport code and applied to 

■ the dynamics of nucleus- nucleus collisions. 

cn 
p , 

1 Introduction 

^ • The many-body theory of strongly interacting particles out of equilibrium is a 

^ I challenging problem since a couple of decades. Many approaches based on the 

• Martin-Schwinger hierarchy of Green's functions [1] have been formulated [2-6] 

and applied to model cases. Nowadays, the dynamical description of strongly 
interacting systems out of equilibrium is dominantly based on on-shell transport 
. , theories and efficient numerical recipies have been set up for the solution of the 

■ coupled channel transport equations [7-13] (and Refs. therein). These transport 
approaches have been derived either from the Kadanoff-Baym equations [14] in 
Refs. [15-19] or from the hierarchy of connected equal-time Green's functions 
[2,20] in Refs. [3,10,21] by applying a Wigner transformation and restricting to 
first order in the derivatives of the phase-space variables {X, p) . 

However, as recognized early in these derivations [3,17], the on-shell quasipar- 
ticle limit, that invoked additionally a reduction of the 87V-dimensional phase- 
space to 7N independent degrees of freedom, where N denotes the number of 
particles in the system, should not be adequate for particles of short life time 
and/or high collision rates. Therefore, transport formulations for quasiparticles 
with dynamical spectral functions have been presented in the past [18] provid- 
ing a formal basis for an extension of the presently applied transport models. 
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Apart from the transport extensions mentioned above there is a further branch 
discussing the effects from non-local collisions terms (cf. Refs. [22-27]). 
The basic questions in this field are: 

i) What are the microscopic mechanisms that lead to thermalization though the 

underlying quantum theory is time reversible? 

ii) What are the effects of off-shell transitions in case of strong coupling? 

iii) Is there an adequate semiclassical limit to the full quantum theory? 

In this article we will attempt to shed some light on these questions for fully 
relativistic self-interacting systems and review some recent progress achieved in 
the last years. Our work is organized as follows: In Section 2 we will study numer- 
ically the self- interacting scalar 0^-theory in 2+1 dimensions on the basis of the 
Kadanoff-Baym equations for homogeneous systems and concentrate on the time 
evolution of correlations and the single-particle spectral function. A comparison 
to the respective Boltzmann limit is performed in Section 3. Furthermore, the 
derivation of a semiclassical off-shell transport theory for inhomogeneous sys- 
tems is presented in Section 4 with special emphasis on the time evolution of 
spectral functions. 

2 The Kadanoff-Baym Equations 

We briefly recall the basic equations for Green's fimctions and particle self- 
energies as well as their symmetry properties that will be exploited throughout 
this work. 

2 . 1 Preliminaries 

Within the framework of the closed-time-path formalism [1,14] Green's functions 
G and self-energies E are given as path ordered quantities. They arc defined on 
the time contour consisting of two branches from (+) -co to oo and ( ) from oo 
to — oo. For convenience these propagators and self-energies are transformed into 
a 2 X 2 matrix representation according to their path structure [4], i.e. according 
to the chronological or antichronological (— ) branch for the time coordinates 
xo and yo- Explicitly the Green's functions are given by 



where the subscript -^y denotes the dependence on the coordinate space vari- 
ables X and y and T"^ (T"") represent the (anti-)timc-ordcring operators. In the 
definition of the factor rj = stands for bosons and 77 = —1 for fermions. 
In the following we will consider a theory for scalar bosons. The full Green's 
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functions are determined via the Dyson-Schwinger equations for path-ordered 
quantities, here given in 2 x 2 matrix representation as 

G> G« J 

/ z'y . 

The self-energies S' are also defined according to their time structure while 
the symbol '©' implies an integration over the intermediate space-time coordi- 
nates from —00 to oo. Linear combinations of diagonal and off-diagonal matrix 
elements give the retarded and advanced Green's functions C^t/adv gg^£_ 

energies S'''^/^^^^ 

(~iret _ /~ic /-^< _ /-^> (~<a /~ia,dv _ /^c /~i> _ fy< (~<a 

xy xy xy xy xy ' xy xy xy xy xy ' 

(3) 

xy xy xy xy xy ' xy xy xy xy xy ' 

Resorting equations (2) one obtains Dyson-Schwinger equations for the retarded 
(advanced) Green's functions (where only the respective self-energies are in- 
volved) 

A— 1 ^ret/adv r , iprct/adv ^ ^rct/adv / a\ 

and the wellknown Kadanoff-Baym equation for the Wightman function , 

^o,l ^xy = -^xz* ^fy + ^xz © ' 

In these equations G~l denotes the (negative) Klein-Gordon differential operator 
which for bosonic field quanta of mass m is given by G~l = —{d^d^ + w?). The 
Klein- Gordon equation is solved by the free propagators Go as 

f Go G<\ /I 0\ 

with the (5-distribution 5xy = 5'''^^^\x — y) for d-t- 1 space-time dimensions. 



2.2 Equilibration Within the Scalar 0^-Theory 

The scalar (/)^-theory is an example for a fully relativistic field theory of interact- 
ing scalar particles that allows to test theoretical approximations without coming 
to the problems of gauge invariant truncation schemes [28-31]. Its Lagrangian 
density is given by (a; = {t, x)) 

C{x) ^ \d,0{x)d"m - A<^4(a;) , (7) 

where m denotes the mass and A is the coupling strength determining the inter- 
action strength of the scalar fields. Its quantization is performed in the canonical 
form (in d spatial dimensions) 

[ct>{x) , dyMU^y, =iS^''\x-y). (8) 
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Since for the real boson theory (7) the relation G^{x,y) = G^{y,x) holds (1) 
the knowledge of the Green's functions G'^{x,y) for all x,y characterizes the 
system completely. Nevertheless, we will give the equations for and G^ 
explicitly since this is the familiar representation for general field theories [5]. 
Self-consistent equations of motion are obtained by a loop expansion of the two- 
particle irreducible (2PI) effective action F on the closed time path. It is given 
by the Legendre transform of the generating functional of the connected Green's 
functions W as 

r[G] = r° + [ln{l-@pGoQpS) + QpGQpE] + $[G] . (9) 

in case of vanishing vacuum expectation value (O|0(x)|O) = [32]. In (9) r° 
depends only on free Green's functions and is treated as a constant while the 
symbols (Dp represent convolution integrals over the closed time path. The func- 
tional <l> is the sum of all closed 2PI diagrams built up by full propagators G; it 
determines the self-energies by functional variation as 

Taking into account contributions up to the 3-loop order the $ functional reads 
explicitly for the problem (7) 

= jj^d^+'x Gi^^xf - ^ 1'^"^'=' l^^^^'y G{x,yf- (11) 

By using the stationarity condition of the action dr/6G = and resolving the 
time structure of the path ordered quantities we obtain the Kadanoff-Baym 
equations for the time evolution of the Wightman function [5,29]: 

- [a^a^+m^] G<{x,y) = E,,^{x) G<{x,y) 

+ jdzojd'^z [S>{x,z)- S<{x,z)\ G<{z,y) (12) 

-Jdzajd^z IJ<{x,z) [G>{z,y)-G<{z,y)], 

- [dp^+m'] G<{x,y) = G<{x,y) 

- jdzojd'^z [G>{x,z)-G<{x,z)\ E<{z,y) 

{x,z) [S>{z,y)-S<{z,y)], 



It 



-Jdzojd'^z 



where d denotes the spatial dimension of the problem (d = 2 in the case consid- 
ered below). Within the 3-loop approximation for the effective action (i.e. the 
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Fig. 1. Self-energies of the Kadanoff-Baym equation: tadpole self-energy (l.h.s.) and 
sunset self-energy (r.h.s.). Since the Unes represent full Green's functions the self- 
energies are self-consistent (see text) 



functional (11)) we get two different self-energies: In leading order of the cou- 
pling constant only the tadpole diagram (l.h.s. of Fig. 1) contributes and leads 
to the generation of an effective mass for the field quanta. This self-energy (in 
coordinate space) is given by 

Sta.d{x) = - iG'^{x,x) , (13) 

and is local in space and time. In next order in the coupling constant (i.e. A^) 
the non-local sunset self-energy (r.h.s. of Fig. 1) enters the time evolution as 

S<'>{x,y) = -^G<'>{x,y)G<'>{x,y)G>'<{y,x) (14) 

— > S<l>{x,y) ^ [G<'>{x,y)]\ (15) 

Thus the Kadanoff-Baym equation (16) in our case includes the influence of a 
mean-field on the particle propagation - generated by the tadpole diagram - as 
well as of scattering processes as inherent in the sunset diagram. 

The Kadanoff-Baym equation describes the full quantum nonequilibrium 
time evolution on the two-point level for a system prepared at an initial time 
ta, i.e. when higher order correlations are discarded. The causal structure of this 
initial value problem is obvious since the time integrations are performed over 
the past up to the actual time xq (or y^, respectively) and do not extend to the 
future. 

In the following we will restrict to homogeneous systems in space. To ob- 
tain a numerical solution of the Kadanoff-Baym equation (12) is transformed to 
momentum space: 

52 G<(p,ti,t2) = -[p'+m2 + Aad(ii)] G<(p,ii,t2) (16) 
- / dt' [S>{pM.t')^S<{pM,t')] G<ip,t',t2) 

Jta 

+ fA'S<{p,ti,t') [G>ip,t\t2)~G<ip,t\h)] 
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where we have summarized both memory integrals into the function 7^. The 

equation of motion in the second time direction t2 is given analogously. In two- 
time and momentum space {p, t, t') representation the self-energies read 

Aad(i) =\ j-^a iG^iP^t^t)^ (17) 
S<{p,t,t') = -- G%,M') G<(r,M') G>ig + r-pXt). 



2.3 Renormalization 

In 2+1 space-time dimensions both self-energies incorporated (1) are ultraviolet 
divergent. Since we consider particles with a finite mass no problems arise from 
the infrared momentum regime. For the renormalization of the divergences we 
only suppose that the time-dependent nonequilibrium distribution functions are 
decreasing for large momenta comparable to the equilibrium ones, i.e exponen- 
tially. Thus we can apply the usual finite temperature renormalization scheme. 
By separating the real time (equilibrium) Green's functions into vacuum (T = 0) 
and thermal parts it becomes apparent that only the pure vacuum contributions 
of the self-energies are divergent. For the linear divergent tadpole diagram we 
introduce a mass counterterm (at renormalized mass m) as 

/(fip 1 
,.^=Vp^+^^ (18) 

that cancels the contribution from the momentum integration of the vacuum part 
of the Green's function. In case of the sunset diagram only the logarithmically 
divergent pure vacuum part requires renormalization, while it remains finite as 
long as at least one temperature line is involved. Contrary to the case of 3+1 
dimensions it is not necessary to employ the involved techniques developed for 
the renormalization of self-consistent theories (in equilibrium) in Refs. [33]. Since 
the divergence only appears (in energy-momentum space) in the real part of the 
Feynman self-energy 1^^ at T = (and equivalently in the real part of the 
retarded/ advanced self-energies iJ^et/adv^^ ^g^j^ absorbed by another mass 
counterterm 

dml^ = -ReS^^oip') = -ReS'^%^''{p') (19) 

f Cpq f (fir 1 U)q + U)r+i^q+r-p 

6 J {2-kY J {2it)'^ AujqWrUJq+r-p [oJq + OJr + OJq+r-p]'^ - pQ 



at given 4-momentum p = {po,p) and renormalized mass m. 
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2.4 Numerical implementation 

For the solution of the Kadanoff-Baym equations we generate a closed set of 
first order differential equations in time for the Green's functions iG^^(x,y) = 
{Hy)(t>{x)), iG<^{x,y) = {(l){y)TT{x)) audiG<^{x,y) = ( 7r(y) 7r(a::) ) with the 
canonical field momentum Tr(x) — dxa4>{x). The disadvantage of having more 
Green's functions in this scheme is compensated by its good accuracy. We es- 
pecially take into account the propagation along the time diagonal which leads 
to an improved numerical precision. The set of differential equations is solved 
by means of a Runge-Kutta algorithm of 4th order. For the calculation of the 
self-energies we apply a Fourier method similar to the one used in [15,34]. The 
self-energies (17) are calculated in coordinate space, where they are products of 
coordinate-space Green's functions ~ that are available by Fourier transforma- 
tion - and are finally transformed to momentum space again. We note that the 
momentum space is discretized by a finite number of momentum states fixed by 
periodic boundary conditions in a finite volume V = a^. 

2.5 Equilibration 

As already observed in the 1-1-1 dimensional case [29-31] the inclusion of the 
mean field does not lead to an equilibration of arbitrary initialized momentum 
distributions, since it only modifies the propagation of the particles by the gener- 
ation of an effective mass. Thermalization in 2-fl dimensions requires the inclu- 
sion of the coUisional self-energies from the sunset diagram. To demonstrate the 




momentum p / m momentum / m 

Fig. 2. Initial Green's functions iC^ {\p\,t = Q,t = 0) (left) and corresponding initial 
distribution functions no{\p\,t — 0) (right) for the distributions Dl, D2 and DT in 
momentum space ( cut of the polar symmetric distribution in px for py = 0) 

phenomenon of equilibration we consider the time evolution of several initial con- 
ditions characterized by the same energy density. The initial equal-time Green's 
functions iG'^{p, t = 0,t = 0) adopted are displayed in Fig. 2 (l.h.s.) as a function 
of the momentum (for py = 0). We concentrate on polar symmetric configu- 
rations due to the large numerical expense for this first investigation^ . Since the 

^ In Section 3 we will present calculations for non-symmetric systems. 



8 W. Gassing and S. Juchem 



0,0 



—1 1 — I — I — I — I I 1 1 



T 1 1 — I — I — r- 



—1 1 — I — I — n- 



distribution D1 
distribution D2 
distribution DT 




p/m=0.8 



p/m=1.6' 



p/m=2.4' 



p/m=3.2 



—1 1 1 1 — I — I I I I 

10 



—1 1 1 — I — I — I I I I 



— I 1 1 — I — I — I I I 



time t m 



100 



1000 



Fig. 3. Time evolution of various momentum modes of the equal-time Green's func- 
tion for \p\/m — 0.0,0.8,1.6,2.4,3.2 (from top to bottom) for three different initial 
configurations Dl, D2, DT (characterized by the different line type) with the same en- 
ergy density. For the rather strong coupling constant \/m — IQ the initial oscillations 
are damped rapidly and have already disappeared at t ■ m = 20. Finally all momentum 
modes assume the same respective equilibrium value independent of the initial state 



equal-time Green's functions (p, t, t) are purely imaginary we display (the 
real part of) i G^. Furthermore, the corresponding initial distribution functions 
no{p,t = 0) are shown in Fig. 2 (r.h.s.). While the initial distributions Dl, D2 
have the shape of (polar symmetric) 'tsunami' waves with maxima at different 
momenta, the initial distribution DT corresponds to a free Bose gas at a given 
initial temperature. The difference between the Green's functions and the distri- 
bution functions is given by the vacuum contribution which has its maximum at 
small momenta, i.e. 2ujpi G^^{p,t~0,t = Q) = 2no(p, t = 0) + l . So even for the 
distributions Dl, D2 the corresponding Green's functions are nonvanishing at 
low momentum. The time evolution of various momentum modes of the equal- 
time Green's function for these three different initial states is shown in Fig. 3. 
Starting from different initial conditions (as shown in Fig. 2) the single momen- 
tum modes finally approach the equilibrium values. For large times all modes 
reach a static limit as characteristic for an equilibrated system. On the other 
hand, for small times one observes a damped oscillating behaviour. This can be 
identified as a typical 'switching-on' effect where the system - being in a static 
(free) situation - is excited by a sudden increase of the coupling constant. In the 
present calculation the initial state is given by a free state with mass m = 1 and 
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the coupling A/m = 16 is switched on for times larger than 0. One might also 

start with an effective initial mass due to the self-consistent tadpole contribution 
[31]. For the study of the equilibration process, however, this does not make any 
significant difference. Due to the character of the self-energies the 'switching-on' 
oscillations are damped in the course of time. The damping depends on the cou- 
pling strength and is much larger for strongly coupled systems. For very strong 
interactions the initial oscillations are even hard to recognize. The final part of 
the time evolution is characterized by an approximately exponential approach 
to the equilibrium value. In contrast to the calculations performed for c/i'^-theory 
in 1+1 space-time dimensions [31] there is no intermediate region, where the 
momentum modes show a power-law behaviour. Furthermore, we observe that - 
depending on the initial conditions and the coupling strength - the momentum 
modes can temporarily exceed their respective equilibrium value. This can be 
seen for the lowest momentum modes of distribution Dl in Fig. 3 which possesses 
initially a maximum at relatively small momentum. Especially the momentum 
mode jpj = of the equal-time Green's function G^, which starts at around 
0.5, is rising to a value of '--^1.8 before it is decreasing again to its equilibrium 
value of about 1.26. Thus the evolution towards the final equilibrium value is - 
after the initial phase with damped oscillations - not necessarily monotonic. On 
the other hand, initial configurations like the distribution DT, where the system 
initially is given by a free gas of particles at a temperature Tq, do not show 
this property. Although they are - of course - not the equilibrium state of the 
interacting theory, they are much closer to it than the distributions Dl and D2. 
Therefore, the evolution towards the equilibrium distribution runs less violently. 



2.6 Time Evolution of the Spectral Function 

Within the Kadanoff-Baym calculations the full quantum character of the two- 
point functions is taken into account. Consequently one fully incorporates the 
spectral properties of the nonequilibrium system during its time evolution. The 
spectral function A{x, y) in our case is given by 

A{x,y) = {[<t>{x) , (l^{y)].) = i [G>ix,y) - G<{x,y)] . (20) 

For each system time T = {ti+t2)/2 the spectral function in Wigner-space is 
obtained via Fourier transformation with respect to the relative time coordinate 
At = h-t2: 

/OO 
dAt e'^^P" A{p,ti=T+At/2,t2 = T-At/2). (21) 
-OO 

A damping of the function A{p,ti,t2) in relative time At corresponds to the 
generation of a finite width of the spectral function in Wigner-space. This width 
in turn can be interpreted as the inverse life time of the scalar particle. We recall 
that the spectral function obeys - for every mean time T and for all momenta 
p — the normalization 

^ Po A{p,po,T) = 1 Vp,T (22) 
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energy p^^ / m energy / m energy p^^ / m 

Fig. 4. Time evolution of the spectral function A{p,pQ,T) for the initial distributions 
Dl, D2 and DT (from left to right) for the two momenta \ p \/m = 0.0 (upper row) and 
\p l/m = 1.6 (lower row). The spectral function is shown for several times t ■ m — 20, 
40, 80, 160, 240, 320, 400 as indicated by the different line types 

which is a reformulation of the equal-time commutation relation (8) . 

In Fig. 4 we present the time evolution of the spectral function for the initial 
distributions Dl, D2 and DT for two different momentum modes \ p\/m = 0.0 
and I p I/to = 1.6. Since the spectral functions are antisymmetric in energy 
A{p, — poi T) — —A{p,po, T) we only show the positive energy part. For our ini- 
tial value problem in two-time space the Fourier transformation (21) is restricted 
for system time T to an interval At G [—ST, 2T]. Thus in the very early phase 
the spectral function assumes a finite width already due to the limited support 
of the Fourier transformation and a Wigner representation is not very meaning- 
ful. We therefore present the spectral functions for various system times starting 
from i • TO = 20 up to i ■ TO = 400. 

For the free thermal initialization DT the evolution of the spectral function is 
very smooth. It is comparable to the smooth evolution of the equal-time Green's 
function for this initialization as discussed in the last Subsection. The spectral 
function is already close to the equilibrium shape at small times being initially 
only slightly broader than for late times. The maximum of the spectral function 
lies (for all momenta) higher than the (bare) on-shcU value and nearly keeps its 
position over the whole time evolution. It results from a positive tadpole mass 
shift which is only partly compensated by a downward shift from the sunset 
diagram. 
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The time evolution of the initial distributions Dl and D2 has a richer struc- 
ture. For the initial distribution Dl the spectral function is broad for small sys- 
tem times (see the line ioi t-m = 20) and is getting a little sharper in the course 
of the system evolution (as presented for the momentum mode \p\/m = 0.0 as 
well as for \ p\/m = 1.6). At the same time the height of the spectral function 
is increasing (as demanded by the normalization property (22)) with time. This 
is indicated by the small arrow close to the peak position. Furthermore, the 
maximum of the spectral function (which is approximately the on-shell energy) 
is shifted slightly upwards for the zero mode and downwards for the mode with 
higher momentum. Although the real part of the (retarded) sunset self-energy 
leads (in general) to a lowering of the effective mass, the on-shell energy of the 
momentum modes is still higher than the one for initial mass m (indicated by the 
'on-sheir arrow) due to the positive mass shift from the tadpole contribution. 
For the initial distribution D2 we find an opposite behaviour. Here the spectral 
function develops at intermediate times a slightly higher width than in the be- 
ginning before it is approaching the narrower static shape. The corresponding 
evolution of the maximum is again indicated by the (bent) arrow. Finally all 
spectral functions assume the (same) equilibrium form. 

As already observed above for the equal-time Green's functions, we emphazise 
that there is no unique type of evolution for the nonequilibrium system. In fact, 
the evolution of the system during the equilibration process is sensitive to the 
initial conditions. 



2.7 Generation of Correlations 

The time evolution within the Kadanoff-Baym equation is characterized by the 
generation of correlations. This can be seen from Fig. 5, where all energy density 
contributions are displayed as a function of time. The kinetic energy density 
£kin is given by all parts of the total energy density that do not depend on 
the coupling constant (oc A°). All terms proportional to are summarized by 
the tadpole energy density etad including the actual tadpole term as well as 
the corresponding tadpole mass coimtcrtcrm. The contributions from the sunset 
diagram (oc A^) - given by the correlation integral as well as by the sunset 
mass counterterm - are represented by the sunset energy density £sun, ^ 

etot(i) =£kin{t) + £tad(i) + esun(i) , (23) 

£tad(t) = I J ^tad(t) iG<^{p,t,t) + i y ^ *nL iG<^{p,t,t) , 

^ Speaking of powers of the coupling constant we mean the 'superficial' order of the 
corresponding diagrams. Since our self-energies E are built up self-consistently by 
full Green's functions G always higher orders of the coupling constant are resummed. 
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Fig. 5. Change of the different contributions to the total energy density in time. The 
sunset energy density £sun decreases as the system correlates. This is mostly compen- 
sated by an increase of the kinetic energy denstty Ekin- Together with the smaller tadpole 
contribution etad the total energy Stot is conserved 



At) 



(2^ 



i/<(p,i,t) 



(2^ 



The calculation in Fig. 5 has been performed for the initial distribution DT 
(which corresponds to a free gas of Bose particles at temperature Tq = 1.59 m) 
with a coupling constant A/m = 16. This state is a stationary solution for the 
well-known Boltzmann equation, but not for the Kadanoff-Baym equation. In 
the latter treatment the system evolves from the uncorrelated state and the cor- 
relation energy density Egun decreases with time. The decrease of the correlation 
energy Esun, which is - with exception of the sunset mass counterterm contribu- 
tion - initially zero, is mainly compensated by an increase of the kinetic energy 
density Ekin- The remaining difference is equal to the change of the tadpole energy 
density etad such that the total energy density is conserved. While the sunset and 
the kinetic term have always the behaviour displayed in the figure, the change 
of the tadpole energy density depends on the initial configuration and can be 
positive or negative. Since the self-energies are obtained within a <P derivable 
scheme (9), the fundamental conservation laws, as e.g. energy conservation, are 
respected to all orders in the coupling constant [32]. 
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3 Comparison to the Boltzmann Limit 

As mentioned above repeatedly the Kadanoff-Baym equation represents the full 
quantum field theoretical dynamics on the one-particle level. However, its nu- 
merical solution is quite involved and it is of strong interest to investigate, in how 
far approximate schemes deviate from the full calculation. Nowadays, transport 
models are widely used in the description of quantum system out of equilib- 
rium (cf. Introduction). Most of these models work in the 'quasiparticlc' picture, 
where all particles obey a fixed energy-momentum relation and the energy is 
no independent degree of freedom anymore; it is determined by the momentum 
and the (effective) mass of the particle. Accordingly, these particles arc treated 
with their ^-function spectral shape as infinitely long living, i.e. stable objects. 
This assumption is rather questionable e.g. for high-energy heavy ion reactions, 
where the particles achieve a large width due to the frequent collisions with 
other particles in the high density and/or high energy regime. Furthermore, this 
is doubtful for particles that are unstable even in the vacuum. The question, in 
how far the quasiparticlc approximation influences the dynamics in comparison 
to the full Kadanoff-Baym calculation, is of general interest [15,34]. 

In order to investigate this question we formulate the Boltzmann limit in 
analogy to Ref . [35] . For the detailed steps and assumptions in the actual deriva- 
tion we refer the reader to Ref. [36] . The Boltzmann equation describes the time 
evolution of the momentum distribution function 

N{p,t) = ^tG<4p,t,t) + ^tG<^{p,t,t) - ReG<^{p,t,t) (24) 
by 2 2 on-shell scattering processes, i.e. 

O^N{p,t) = /-^ (27r)'^,5W(p+q-r-s) (25) 

{NpNqNrNs - NpNqNrNs} IT 5{LJp + UJq-LJr -i^s) 

using Np = N{p, t) and Np = N{p, t) + 1 for the corresponding Bose factors. 

The numerical procedure for the solution of (25) is basically the same as 
the one developed for the solution of the Kadanoff-Baym equation. Moreover, 
we calculate the actual momentum dependent on-shell energy iVp for every mo- 
mentum mode by a solution of the dispersion relation including contributions 
from the tadpole and the real part of the (retarded) sunset self-energy. Thus it 
is guaranteed, that the particles are treated as quasiparticles with the correct 
energy-momentum relation at every time. 

For the comparison between the full Kadanoff-Baym dynamics and the Boltz- 
mann approximation we concentrate on equilibration times. As a measure for 
the equilibration time we consider the time scale on which the initially non- 
isotropic distribution proceeds to the polar symmetric equilibrium value. We 
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Fig. 6. Relaxation rate (divided by the coupling squared) for Kadanoff-Baym and Boltz- 
mann calculations as a function of the interaction strength. For the two different initial 
configurations the full Kadanoff-Baym evolution leads to a faster equilibration 



define a moment Q{t) for a given momentum distribution n{p, t) at time t by 



(27r)2_U^2__fV 



= \- a^, : ' (26) 



which vanishes for symmetric systems, e.g. for the equiUbrium state. For the 
Kadanoff-Baym case we calculate the actual distribution function by 

n{v,t) = ^G<^{p,t,t) G<^{p,t,t) - i . (27) 

The moment Q{t) shows an approximately exponential decrease in time such 
that we can define a relaxation rate Fq via the relation Q{T) oc e~^'^*. In Fig. 
6 the relaxation rate Fq (scaled by the coupling strength squared) is displayed 
for the Kadanoff-Baym and the Boltzmann calculation for two different initial 
configurations. The initial distribution 2 corresponds to the initial state of two 
(on the Paj-axis) separated particle accummulations (cf. Fig. 7). The peak at low 
momenta again comes from the vacuum contribution of the Green's function. The 
time evolution within the full Kadanoff-Baym theory is presented in Fig. 7 by 
several snapshots at times t-m = 0, 10.5, 21, 35, 70, 105 showing the propagation 
towards an isotropic (static) final state. For distribution 1 the position and the 
width of the two bumps is modified. 
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Fig. 7. Evolution of the Green's function in momentum space. The equal time Green's 
function is displayed for various times t ■ m — 0, 10.5, 21, 35, 70, 105. Starting from 
an initially non-isotropic shape it develops towards a symmetric final distribution 

Fig. 6 shows for both initializations that the relaxation in the full quantum 
calculation happens faster for large coupling constants than in the quasi-classical 
approximation, whereas for small couplings the equilibration times of the full and 
the approximate evolutions are comparable. While the scaled relaxation rate 
Fq/X^ is nearly constant in the Boltzmann case, it increases with the coupling 
strength in the Kadanoff-Baym calculation (especially for initial distribution 2). 

Since the free Green's function - as used in the Boltzmann calculation - has 
only support on the mass shell, only (2 ^ 2) scattering processes are described 
in the Boltzmann limit. All other processes with a different number of incoming 
and outgoing particles vanish. Within the full Kadanoff-Baym calculation this 
is very much different since here the spectral function ~ determined from the 
self-consistent Green's function - aquires a finite width. The Green's function 
has support at all energies - although it drops far off the mass shell. Especially 
for large coupling constants, where the spectral function is sufficiently broad, the 
three particle production process gives a significant contribution to the collision 
integral. Since the width of the spectral function increases with the interaction 
strength, such processes become more important in the high coupling regime. 
As a consequence the difference between both approaches is larger for stronger 
interactions as observed in Fig. 6. For small couplings X/m in both approaches 
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basically the usual 2 <-> 2 scattering contributes and the results for the thermal- 
ization rate Fq arc quite similar. 

In summarizing this Section we point out that the full solution of the Kadanoff- 
Baym equations does include 4, 1 3 and 2 2 off-shell collision processes 
which - in comparison to the Bolzmann on-shcU 2^2 collision limit ~ become 
important when the spectral width of the particles reaches about 1/3 of the 
particle mass. On the other hand, the simple Boltzmann limit works surpris- 
ingly well for smaller couplings and those cases, where the spectral function is 
sufficiently narrow. 

4 Derivation of Semiclassical Transport Equations for 
Particles with Dynamical Life Times 

The Kadanoff-Baym equations (5) presently cannot be solved for general inho- 
mogeneous problems due to the high complexity of these equations and the huge 
'computer storage' required. It is thus of interest to obtain a semiclassical limit 
which is much easier to solve e.g. by current testparticle methods. 

We start again with the Kadanoff-Baym equation (5), however, change to the 
Wigner representation via Fourier transformation of the rapidly oscillating rela- 
tive coordinate {x — y). The theory is then formulated in terms of the coordinates 
X = {x + y)/2 and the momentum p, 

Fxp = J d\x-y) e^fMC-^-f") F^y. (28) 

Since convolution integrals convert under Wigner transformations as 

j d\x-y) e^PM(-^-yn F,,,,QF2,,y = e'^^ i^2,Xp, (29) 

one has to deal with an infinite scries in the differential operator o which is a 
four-dimensional generalization of the Poisson- bracket, 

As a standard approximation of kinetic theory only contributions up to first 
order in the gradients gradient expansion are considered. This is justified if the 
gradients in the mean spatial coordinate X and momentum p are small. 
Applying this approximation scheme (Wigner transformation and neglecting all 
gradient terms of order n > 2) to the Dyson-Schwinger equations of the retarded 
and advanced Green's functions one ends up with 

(/ - m2 - ReSTp )ReGT^ = l-l Fxp Axp , 

( p2 - - ReSTp ) Axp = Fxp Re G% , (31) 
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where we have separated the retarded and advanced Green's functions as well 
as the self-energies into real and imaginary contributions 

The imaginary part of the retarded propagator is given (up to a factor) by 
the normalized spectral function Axp (20) while the (negative) imaginary part 
of the self-energy is half the width Fxp- Above equations (31) arc given for 
a (renormalized) vacuum mass Mq and a (renormalized) retarded self-energy 
S'^^*. In the presence of an additional space-time local self-energy (e.g. tadpole 
self-energy) the mass term is shifted accordingly Mq -|- Stad- From the 

algebraic equations (31) we obtain a direct relation between the real and the 
imaginary part of the propagator (provided Fxp 7^ 0): 



ReGTp = - ^Axp. (33) 

' Xp 



The algebraic solution for the spectral function shows a Lorentzian shape with 
space-time and four- momentum dependent width Fxp- This result is valid for 
bosons to first order in the gradient expansion, 



Fxp 

{p^-M^-ReS%y + Fy4 



For the real part of the retarded Green's function we get also algebraically 

ret - ^0 - R-eZ-f t 

^ ^xp = 7-2 J72 5rF^in2~~r^^2~~77 • 



4.1 Transport Equations 

The Kadanoff-Baym equation (5) gives in the same semiclassical approximation 
scheme a generalized transport equation, 

o{p^-MS-ReErp}{G<^} - o{S<^} {ReG%} 

= 2 [ ^xp Gxp - ^xp G^p ] > (36) 

and a generalized mass-shell equation, 

[p2 - Ml - R^ETp] G% - S<^ ReG% 

= \o{S<^}{Axp} - \o{Fxp}{G<x^}. (37) 

In the transport equation (36) one recognizes on the l.h.s. the drift term 'p^'d^* 
as generated by the contribution o{p^ — Mo}{»}, as well as the Vlasov term 
determined by the real part of the retarded self-energy. On the other hand the 
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r.h.s. represents the collision term with its 'gain and loss' structure. To evaluate 

the o{Z'<}{Rc G"'''*}-tcrm in (36), which docs not contribute in the quasiparticlc 
limit, it is useful to introduce distribution functions for the Green's functions 
and self-energies as 

i G<p = Nxp Axp, * Gxp = ( 1 + ^xp ) Axp , 

i^xp = N^pTxp, iS^p = (1 + N^p)rxp. (38) 

Following the argumentation of Botcrmans and Malflict [16] the distribution 
functions N and in (38) should be equal in the second term of the l.h.s. 
of (36) within a consistent first order gradient expansion. As a consequence the 
self-cncrgy E< can be replaced by ■ F/A in the term «{i7<}{ReG''''*}. The 
general transport equation (36) then can be written as 

o{p'-M^-EeSTp}{G<^} - -^o{rxp}{{p'-MS-ReST,)G<„} 

X AxpTxp = i [S^pG<^ - S^pG>^] . (39) 



4.2 Test Particle Representation 

In order to obtain an approximate solution to the transport equation (39) we 
use a test particle ansatz for the Green's function more specifically for the 
real and positive semidefinite quantity 

N 

Fxp = AxpNxp =^G<^^Y. S^^\X-X,{t)) 6^^\p-pS)) 6{po-e^{t)). (40) 

In the most general case (where the self-energies depend on four-momentum 
p, time t and the spatial coordinates X) the equations of motion for the test 
particles read 



dt 

dpi 
dt 

dcj 
dt 



1 — G(i) 2ej 



-1 1 



2pi + Vp,ReSJ^^ + 



el-p^-M§-ReS^^^ 



1 ~ G{i) 2ej 



1 - C(i) 2e 



Vx.Rerr* + ^— ^ ^ VxT^i) 



dt 



dt 



(41) 
(42) 
(43) 



where the notation implies that the function is taken at the coordinates of 
the test particle at time t, i.e. = F{t, Xi{t),Pi{t),ei{t)). 
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In (41-43) a common multiplication factor (1 — C(i)) ^ appears, which con- 
tains the energy derivatives of the retarded self-energy 



(44) 



It yields a shift of the system time t to the 'eigentime' of particle i defined by 

ti = t/{l — C(i))- As the reader immediately verifies, the derivatives with respect 
to the 'eigentime', i.e. dXi/dti, dpi/dti and dci/dti then emerge without this 
renormalization factor for each test particle i when neglecting higher order time 
derivatives in line with the semiclassical approximation scheme. In the limiting 
case of particles with vanishing gradients of the width Fxp these equations of 
motion reduce to the well-known transport equations of the quasiparticle picture. 

Following Refs. [37] we take = — Re 17''''* as an independent variable 
instead of po) which then fixes the energy (for given p and M2) to 

pI= p-" + + Re 1:5,%.. (45) 

Eq. (43) then turns to 

d(Mf - M2) _ Mf- Mi dr^^ 



dt r'fi) dt 



(46) 



for the time evolution of the test particle i in the invariant mass squared as 
derived in Refs. [37]. For applications of the semiclassical off-shell transport 
approach we refer the reader to Refs. [37]. 



5 Summary 

In this article we have studied the time evolution of an interacting field theo- 
retical system, i.e. (/)'*-field theory in 2+1 space-time dimensions, on the basis 
of the Kadanoff-Baym equations for a spatially homogeneous system including 
the self-consistent tadpole and sunset self-energies. We find that equilibration 
is achieved only by inclusion of the sunset self-energy. Simultaneously, the time 
evolution of the single-particle spectral function has been calculated for various 
initial conditions. A comparison of the full solution of the Kadanoff-Baym equa- 
tions with the solution for the corresponding Boltzmann equation shows that a 
consistent inclusion of the spectral function has a sizeable impact on the equili- 
bration rates if the width of the spectral function becomes larger than ~l/3 of 
the particle mass. 

Furthermore, the conventional transport of particles in the on-shell quasipar- 
ticle limit has been extended to particles of finite life time by means of a dy- 
namical spectral function A{X,p,M'^). Starting again from the Kadanoff-Baym 
equations we have derived in consistent first order gradient expansion equations 
of motion for test particles with respect to their time evolution in X,p and 
M^. This off-shell propagation has been examined for a couple of model cases in 
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Refs. [37] as well as for nucleus-nucleus collisions showing that - at subthreshold 
energies - the off-shell dynamics play an important role for the production of 
energetic particles. 
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